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Measurements have been made of the probability distribution of total transmission of 
microwave radiation in waveguides filled with randomly positioned scatterers which would 
have values of the dimensionless conductance g near unity. The distributions are markedly 
non-Gaussian and have exponential tails. The measured distributions are accurately de- 
scribed by diagrammatic and random matrix calculations carried out for nonabsorbing 
samples in the limit <? S> 1 when g is expressed in terms of the variance of the distribution, 
which equals the degree of long-range intensity correlation across the output face of the 
sample. 
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Nonlocal correlation in the flux transmitted through mesoscopic samples leads to enhanced fluctuations 
of local and spatially averaged transmission for both classical and quantum waves. Such fluctuations 
increase dramatically as the ensemble average of the dimensionless conductance, g, approaches unity. Low 
values of g can be achieved in quasi-one-dimensional samples such as conducting wires or multimode waveg- 
uides with lengths much greater than the transverse dimensions. In this Letter, we report measurements 
of the probability distribution of total transmission of microwave radiation in long waveguides filled with 
randomly positioned scatterers which in the absence of absorption would have values of g near unity. 
The distributions observed are markedly non-Gaussian. They are compared to recent diagrammatic and 
random- matrix calculations for nonabsorbing samples in the limit g ^ 1 . This is done by reexpressing 
the distribution, which is a function of the single parameter g, as a function of the variance of the normal- 
ized transmission using the relation between these parameters. This result is in excellent agreement with 
the measured transmission distributions and indicates that the variance of the normalized transmission, 
which equals the degree of long-range intensity correlation across the output face of the sample, is the 
essential parameter describing fluctuations in random media. 

Key transmission quantities in order of increasing spatial averaging are the intensity. Tab, which is the 
transmission coefficient for incoming mode a into mode 6, the total transmission for incoming mode a, 
Ta = X^b^ab ~ ^/L, and the total transmittance T = X^ah-^ab ^ Nl/L, where I is the transport mean 
free path, L is the sample length, and N is the number of modes. The total transmittance is equivalent 
to the dimensionless conductance in electronic systems, T — Gjie^jK), where G is the conductance |^ and 
g = < T > — Ni/L. Though the variances of the transmission quantities normalized to their ensemble 
average values, Sab — Tab/<Tab>, Sa — Ta/<Ta> and s = T/<T>, are reduced as the extent of spatial 
averaging increases, fluctuations in these quantities do not self average, as they would if spatial correlation 
were absent. To leading order in l/g, the enhancement of the variances of Sab, Sa, and s arising from 
nonlocal correlation is 1, L/£ and {L/lY respectively, which results in values of the variances of 
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To examine the scaling and the universality of transport, it is important to measure the full distribution 
of key transmission quantities as the sample size, and hence g, changes. In previous work, nonlocal 
correlation has been shown to lead to higher probabilities at large values of the intensity, leading to a 
deviation from negative exponential statistics for polarized microwave radiation when g ^ \Q | pd] , p^ , as 
well as to discernable deviations from a Gaussian distribution and enhanced variance for the total optical 
transmission when g > 10"^ [Q. 

Recently, an expression for P{sa) in terms of g for nonabsorbing samples was obtained by Nieuwenhuizen 
and van Rossum using diagrammatic techniques combined with random matrix theory ^ and subsequently 
by Kogan and Kaveh within the framework of random matrix theory Q. The diagrammatic calculations 
neglect some terms of order higher than l/g, whereas computations based on random matrix theory neglect 
sample-to-sample fluctuations in the probability distributions of eigenvalues of the transmission matrix and 
are expected to be accurate only to order 1/g. More recently, van Langen, Brouwer and Beenakker carried 
out a nonperturbative calculation of the total transmission distribution in the absence of absorption . 
An analytic solution is obtained for the case in which time reversal invariance is broken {[3 — 2) but not for 
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the case of time reversal symmetry (/? = 1) considered here. However, good agreement is found between 
the /3-independent result for P{sa) obtained in Refs. and the result for /? = 2 in Ref. for g > 10. 

The distribution of total transmission has been measured previously by de Boer et al. in optical mea- 
surements in slabs of titania particles jl^ . Samples with g > 10^ were studied and the distribution was 
found to be Gaussian to within 1%. A measure of the deviations of the distribution from a Gaussian is the 
value of the third cumulant <sl>c which gives the skewness of the distribution and vanishes for a Gaussian 
distribution. For the samples studied, <s^>c was of order of 10~^. It was found that <sj^>c= 7g <Sa>c 
with 7g = 2.9 ± 0.4 which is consistent with the value calculated for Gaussian beam excitation of 3.20 

In the present work, low values of g are achieved by placing the sample in a cylindrical copper tube in 
order to restrict transverse diffusion and thus the number of modes N. The samples consist of randomly 
positioned polystyrene spheres with diameters of 1.27 cm at a volume filling fraction / = 0.55. Transmission 
spectra were taken at tube diameters of 7.5 and 5.0 cm and various sample lengths in the frequency range 
16.8 - 17.8 GHz. The microwave radiation is coupled to the sample by a 0.4 cm wire antenna placed 
0.5 cm from the front surface of the sample. The frequency is incremented in 4 MHz steps. The sample 
tube is rotated between successive measurements to produce new scatterer configurations. The total 
transmission is measured by use of a single Schottky diode detector positioned inside an integrating sphere 
which rotates at 2 Hz around the sample axis. The integrating sphere has a diameter of 40 cm and is 
comprised of two concentric plastic spherical shells separated by a distance of 2 cm. The outer shell is 
covered with aluminum foilt to form an irregular reflecting surface. The region between the shells is filled 
with thin-walled aluminum cylinders with diameters of 0.75 cm and typical lengths of 1 cm. The cylinders 
tumble as the integrating sphere rotates, resulting in fluctuations of the intensity at the detector with a 
correlation time of 2 ms for the sample with a length of 100 cm. The signal is averaged for 1 s at each 
frequency, giving an uncertainty of 2.5 % in the measurement of transmission. The signal is normalized 
by its ensemble average to give Sa- The transmission distributions P{sa) are obtained by using the data 
from at least 1000 sample configurations. Distributions obtained using different intervals of the frequency 
range coincide within experimental error. In the frequency range of the measurements, £ 5 cm and 
N = k'^(P/8 « 200 and 90 for samples in tubes with diameters d — 7.5 and 5.0 cm, respectively. The wave 
number k = 2tt/X ~ 2T:vn/c is calculated using an effective medium index of refraction n « 1.4. 

The transmission distributions for three samples with dimensions (a) d — 7.5 cm, L — 66.7 cm, (b) 
d ~ 5.0 cm, L = 50.0 cm, and (c) d = 5.0 cm, L — 200 cm are shown in Fig. 1. In the absence of 
absorption, the dimensionless conductance for these samples without localization corrections, g — N£/L, 
would be approximately 15.0, 9.0, and 2.25 for samples a, 6, and c, respectively. The distribution broadens 
and the deviation from a Gaussian becomes more pronounced as either the sample length increases or 
the cross-sectional area decreases. A value of <s^>c as large as 0.112 ± 0.003 is observed for sample c. 
Deviations from a Gaussian distribution in the tail of the distribution for this sample can be seen in the 
semilog plot of P{sa) in Fig. 2. For values of > 2 the distribution is nearly exponential. 

In Fig. 3, we present a plot of <s^>c versus <s^>^. The solid line is a least square linear fit to the data 
which gives 7 = 2.38 ± 0.05. Within experimental error this equals the value 7 — 2.40 calculated for an 
incident plane wave in the lowest order of a diagrammatic perturbation expansion in the small parameter 
1/g 1^. The results are compared to calculations for a plane wave since d <^ L and there is a nearly 
complete mixing of modes in the sample, giving a uniform average intensity along a cross section of the 
sample. The agreement between theory and experiment is surprising, however, since l/g ^ 0.1 for all 
samples, reaching a value of approximately 0.3 for sample c, and is by no means small. Furthermore, the 
influence of absorption was not included in the calculations, whereas the samples used in the experiment 
are strongly absorbing with L > « 40 cm, where La is the exponential absorption length |l5j . 

We now consider the full transmission distribution. The theoretical expressions for the full distribution 
function in Refs. are given as functions of g for nonabsorbing samples. In the present case of 

strong absorption, the photon number is not conserved, and g cannot be defined in terms of the steady 
state transmission, while serving as a useful measure of the proximity to the localization transition. This 
can be seen by noting that the reduction of the average transmission due to absorption would lead to a 
reduced value of g even though the presence of absorption tends to lessen the degree of correlation in the 
sample and to push the system farther from the localization threshold. On the other hand, a parameter 
which characterizes the transmission distribution as well as the closeness to the localization threshold, 
even in the presence of absorption, is the degree of correlation of intensity in different coherence areas of 
the transmitted speckle pattern, <SsabSsab' >■ Were this correlation to vanish, fluctuations in different 
coherence areas would be independent and the transmission distribution would be Gaussian by the central 
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limit theorem with var{sa) = <s^>c= l/N. As a result of nonlocal correlation, however, the variance of 
the transmission is enhanced. It is given by <s^>c= (<s^j>c — 1)/2 = <SsabSsab'> |3,P, [To| , p^ ■ The 
last equality is consistent with the results of Ref. JlGf when the cumulant intensity correlation function is 
properly normalized to the renormalized average transmission [0. In that case, the crossing parameter 
found by Shnerb and Kavch [ p8[ which determines the intensity distribution is found experimentally to be 
equal to <SsabSsab'> The connection of <6sabSsab'> to the full transmission distribution can be 

seen by considering the expression of Rcfs. |^,^ for P{sa) in the absence of absorption in the limit g ^ 1, 

/■+'°° dx 

PiSa)^ —exp[xSa - ^{x)], (1) 

J-ico 27r« 

where 

= gln'^ {y/l + x/g + \/x/g) 
is the generating function. From Eq. (|l|) one obtains the expression for <s^>c in terms of g, 

<4>c^ |-- (2) 

From these expressions, a general relation for P{sa) in terms of <s^>c, or equivalently <SsabSsab'>, can 
be found by using Eq. (||) to define a new parameter g' = 2/3<s^>c which is substituted for g into Eq. (|^). 
Plots of P{sa) obtained by following this procedure with g' determined from the measured values of <s^>c 
are shown as the solid lines in Figs. 1 and 2. We find that P{sa) is accurately given even for the lowest 
value of g' of 3.06 (sample c). The distribution of Eq. (|l|) with g' substituted for g gives the exponential 
tail, P{sa) ^ &xp{—g' So) in the limit Sa ^ 1. For Sa > 2.0, the linear fit to the logarithm of the measured 
transmission distribution for sample c gives a slope of 2.71 ± 0.06 in accord with the exponential fit of 
the theoretical curve of 2.70 in this range and is quite close to its predicted asymptotic value of 3.06 for 

Sa > 1- 

The extent of the agreement of Eq. (0) when g' is substituted for g can also be gauged from the comparison 
between the calculated (circles) and the measured (squares) moments of the transmission distribution shown 
in Fig. 4 for samples with g' — 10.2 ±0.1 and g' = 3.06 ± 0.04. The moments calculated from the theory 
are close to those obtained from the measured distributions. At n — 10, these defer by approximately 
10% which is within the experimental error. Thus it appears that P{sa) can be expressed as a function of 
the parameter <s^>c. 

The agreement between theory and experiment indicates that the ratio of moments is accurately reflected 
in Eq. (|l|). The dependence of the second cumulant itself upon sample dimensions is shown in Fig. 5. In 
the limit g ^ 1, in the absence of absorption, <s^>c= 2L/3NI. The straight line in the figure is drawn 
through the first data point and represents <s^>c~ L/N. As g — > 1, and the localization threshold 
is approached, the scaling theory of locahzation [|l9| suggests that g falls more rapidly and hence <Sq>c 
should increase superlinearly with sample length. Instead, we find that <s^>c depends sublinearly upon L. 
This is presumably due to the presence of absorption which diminishes the degree of nonlocal correlation. 
This raises the question of whether the transmission distribution continues to broaden as L increases or, 
instead, it reaches a limiting distribution for particular sample parameters. 

In conclusion, we have measured the total transmission distribution of microwave radiation in quasi- 
one-dimensional absorbing samples with small values of g. We find that the distribution can be described 
using an expression originally derived for nonabsorbing samples in the limit g ^ 1 when this expression is 
reformulated as a function of the single parameter g' = 2/3 <s^>c determined from the measurements. The 
validity of the expression for values of g' as small as 3, well beyond the limits assumed in the calculations, 
may well be associated with the identification of <Sa>c with <Ssab6sab'>, the degree of spatial correlation 
in the sample, which is the key microscopic parameter in mesoscopic physics. 
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FIGURES: 



Fig. 1. Distribution function of the normalized transmission P{sa) for three samples with dimensions (a) 
d = 7.5 cm, L = 66.7 cm, (b) d = 5.0 cm, L = 50.0 cm, and (c) d = 5.0 cm, L = 200 cm. 

Fig. 2. Semilogarithmic plot of the transmission distributions for the same samples as in Fig. 1. 

Fig. 3. Plot of <s^>c versus <s^>^. The solid line represents a least square linear fit to the data. 

Fig. 4. Comparison of the calculated (circles) and measured (squares) moments of the transmission distri- 
bution for samples with (a) g' = 10.2 and (b) g' = 3.06. 

Fig. 5. Dependence of the second cumulant <s„>c upon sample dimensions. 
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(a) g' = 10.2 
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(b) g' = 3.06 
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